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Abstract We present a set of restrictions on the fractional differential equa- 
tion x^(t) = g(x(t)), t > 0, where a G (0,1) and g(0) = 0, that leads to 
the existence of an infinity of solutions starting from x(0) = 0. The operator 
x^ is the Caputo differential operator. 
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1 Introduction 



The issue of multiplicity for solutions of an initial value problem that is 
associated to some nonlinear differential equation is essential in the modeling 
of complex phenomena. 

Typically, when the nonlinearity of an equation is not of Lipschitz type 
PP, there are only a few techniques to help us decide whether an initial value 
problem has more than one solution. As an example, the equation x' = 
f(x) = y/x- X(o,+oc)(x) has an infinity of solutions x T (t) = (t ~ 4 T) ■ X(t,+oo)(*) 
defined on the nonnegative half-line which start from x(0) = 0. Here, by \ 
we denote the characteristic function of a Lebesgue- measurable set. 

An interesting classical result j9j[T], which generalizes the example, asserts 
that the initial value problem 



where the continuous function g : R — > R has a zero at x and is positive 
everywhere else, possesses an infinity of solutions if and only if J xq+ < 
+00. 

Recently, variants of this result have been employed in establishing various 
facts regarding some mathematical models [121 02]. In particular, if the 
function g is allowed to have two zeros x < x± while remaining positive 
everywhere else and 



then the problem ([T]) has an infinity of solutions (xt)t>o such that lim xr{t) 



Our intention in the following is to discuss a particular case of the above 
non-uniqueness theorem in the framework of fractional differential equations. 
To the best of our knowledge, the result has not been established in its full 
generality yet for any generalized differential equation. We mention at this 
point the closely connected investigation [I]. 

In the last number of years, it became evident that differential equations 
of non-integer order, also called fractionals (FDE's), can capture better in 
models many of the relevant features of complex phenomena from engineer- 
ing, physics or chemistry, see the references in [21, [3J [5j [HJ HH [151 QUI El E] . 



Let us consider a function h G C X {I, R) n C(I, R) with limjt 1 - /^)] G R 
for some a G (0, 1), where / = (0, +00). The Caputo derivative of order a of 




x'(t) =g(x(t)), t>0, 

x(0) = Xq, Xq G R 



(1) 




t— >+oc 



= Xi. 
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h is defined as 



r(l-a) J (t 

where T is Euler's function Gamma, cf. [T5J p. 79]. To have an idea about the 
eventual smoothness of this quantity, we mention that by letting the function 
h' be at least absolutely continuous [161 Chapter 7] the (usual) derivative of 
will exist almost everywhere with respect to the Lebesgue measure m 
on M, see [T7] p. 35, Lemma 2.2]. Further, we have that 

provided that hS a ' is in L°°(m). 

The initial value problem we investigate in this paper is 



x^(t)=g(x(t)), tel, 
x(0) = 0, 



(3) 



where the function g : R — > R is continuous, g(0) = and g(u) > when 
u G (0, 1]. Further restrictions will be imposed on g to ensure that f Q+ < 
+oo. 

By means of (j2J), we deduce that 
= /" x'(s)ds 

dr ds 



o 



r(a)7 Lr(l-a)J (a-r 

and so the problem ([3]) can be recast as 



J V(t) = 9 (H , t > 0, 

I y(0) = o, 



(4) 



where y = x^ a \ (3 = 1 — a and the (general) function g has absorbed the 



constant ^7-^. 

r(a) 



In the next section, we look for a family (ut)t>o, with yx G C([0, 1],' 
of (non-trivial) solutions to (j3J). 
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2 Infinitely many solutions to (ED) 

We start by noticing that the function / : [0, 1] — > (0, 1) with the formula 
f( x ) = 2+X-/3 = ^ ~ 2+x-p * s increasing. Introduce now the numbers 8\ and 
5 2 such that 

/(l) = 3^ > ft > ft > /(0) = (5) 

Obviously, Si, S 2 G (0, 1). 

Set ci < C2 in (0, +oo) and assume that 

ci-y h <g{y)<c 2 -y & \ yE [0,1]. (6) 

As a by-product, < -f^^ < +oo. Further, suppose that there exists 

c > such that 

\g{y2) -g{yi)\ < —— ■ \vi -yi|, 2/1,2/2^(0,1]. (7) 

(min{j/i,2/2» 

The latter condition has been inspired by the analysis in [14] . 
Introduce now the numbers Yi,Y 2 >_\ and T G (0, 1) with 

(Y 1 + Y 2 )(l -Tf-P < 1 -/? (8) 

and 

8Yi < Cl < c 2 < (1 - P)Yt 52 (9) 

and 

vt) <l (10) 

These will be used in the following for describing the solution y?- 

Several simple estimates, of much help in the proof of our result, are 
established next. Notice first that, via the change of variables s = T+u(t—T), 
we get 



*(a-T) 



l+e r 1 7|!+ e 



= J B(2 + e,l-/3)(t-T) 2+£ - /3 , (11) 
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where e G (0, 1) and B represents Euler's function Beta [T5]. Also, 

" J du 



B(2 + e,l-fi) > 



1 u 1+e , 1 
du > 



l (1 -u)P ~ 2!+ £ Ji (l-uf 

2 -(2+e-/3) 1 
> 



1-/3 -8(1-/3) 



and 



s(2+ £ ,i-/3)< r 

./Q 



(1-m)/ 3 1-/3' 

Now, returning to (jHJ), there exist S\, e 2 G (0, 1), with £1 > £ 2 , such 

f(e 1 ) = 5 1 >5 2 = f(e 2 ). 

In particular, 1 — 81 = 2 +7f-p anc ^' means of (|T3l) . 

KiS(2 + £i,l- /3)(1 - T) 2+£l "^ 
+Y 2 5(2 + £ 2 , 1 - /3)(1 - T) 2+£2_/3 < 1. 

Taking into account f[T2"j) and (|UJ), we deduce that 

Si 



which leads to 



|8(l-0) J 
Yi < ci[y 1 S(2 + e 1 ,l 



< Cl 5(2 + £!,!-/?) 



<5l 



1* 



and, via (113]) . 



c 2J B(2 + £ 2 , 1 - /3) 52 < . ° 2 .. < ' " 



" 1-/3' 
which implies that 

C2 [y 2 £?(2 + £ 2 ,i-/3)] 52 <y 2 . 

Let the set y C C( [T, 1] , M) be given by the double inequality 

Y x {t - T) 1+£l < y(t) < Y 2 (t - T) 1+£2 , * G [T, 1], y G y. 
Observe that, by means of (TIB]). ffTTj) . (JB}, 

" y(a) 



7^(2 + £1, 1 - /3)(£ - Ty +£l -v < 



(t - 8 )P 



ds 



and 

Yi(t-T) 1+£l < c 1 [F lJ B(2 + £ 1 ,l-/3)]' 51 • (t-T) 1+£1 

= Cl [Y 1 B(2 + e 1 ,l-m-T) 2+£ ^} 51 (20) 

< ,(T A* 



(t - s y 

Similarly, via (IT7|) . 
and 

^(/(t^ ds ) ^ c 2 [Y 2 B(2 + e 2 ,l-(3)] 5 ^(t-T) 1+ ^ 
< Y 2 (t-T) 1+£2 . 

In conclusion, the mapping t\-^gyj^, (t-s)P * s a mem ber of 3^ when- 
ever y E y. Also, taking into account (Tl5]) . we deduce that the quantities y = 
yiS(2+ei, l-/?)(t-T) 2+ei -0 from (T50D and y = Y 2 B(2+e 2 , l-/3){t-T) 2+e2 ^ 
from fT2~T]) belong to [0, 1] as imposed in (jSJ). 

We are now ready to state and prove our main result. 

Theorem 1 Assume that the nonlinearity g of (fJJ) satisfies the restrictions 
(EP; (0)- Given the numbers Y x , Y 2 , T subject to (TJ|) ; / TIflj) and the 
set y = y(Yi,Y 2 ,T) from < T73)) . the problem (T^J) /ias a unique solution y? in 

y. 

Proof. The operator O : y — > y with the formula 

0(y)(t) = g (£ j0^dsY yey,te[T,l], 
is well defined. 

The typical sup-metric d(yi, y 2 ) = sup \yi(t) — y 2 (t) \ provides the set 3^ 

te[T,i] 

with the structure of a complete metric space. 
Taking into account (J7|), (|T9|) and ({TO]) , we get 

|0(lfc)(t)-0(lfc)(*)| 

< 



[Y l B(2 + e 1 ,l-f3)(t-T)^-P] 1 - 51 Jt (t - s)? 



< c 



k-d(y!,y 2 ), yi, y 2 G 3^- 
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The operator O being thus a contraction, its fixed point i/t in y is the 
solution we are looking for. Notice that x/t is identically null in [0,T]. □ 
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